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Abstract 



The mass of a black hole is interpreted, in terms of thermodynamic potentials, 
as being the enthalpy, with the pressure given by the cosmo logical constant. 
The volume is then defined as being the Legendre transform of the pressure 
and the resulting relation between volume and pressure is explored in the case 
of positive pressure. A virial expansion is developed and a van der Waals 
like critical point determined. The first law of black hole thermodynamics 
includes a PdV term which modifies the maximal efficiency of a Penrose 
process. It is shown that, in four dimensional space-time with a negative 
cosmological constant, an extremal charged rotating black hole can have an 
efficiency of up to 75%, while for an electrically neutral rotating back hole 
this figure is reduced to 52%, compared to the corresponding values of 50% 
and 29% respectively when the cosmological constant is zero. 



PACS nos: 04.60.-m; 04.70.Dy 



1 Introduction 



The thermodynamics of black holes is a rich and fascinating area of research 
which continues to yield surprises. The first law of black hole thermodynam- 
ics is usually written as 

dM = TdS + VLdJ + $dQ (1) 

where T = ^- is the Hawking temperature of the black hole (with k the 
surface gravity), S = j the entropy (with A the area in Planck units), f2 
the angular velocity, J the angular momentum, $ the electrostatic potential 
difference between infinity and the horizon, Q the electric charge and M 
the mass. The mass is usually interpreted as the internal energy, in the 
thermodynamic sense, of the black hole, but it was suggested in [I] that it 
is more correctly interpreted as the enthalpy. In this context it is notable 
that there is no PdV term in (CQ), corresponding to a change in volume at 
ambient pressure P. When a cosmological constant, A, is included there is 
a natural candidate for a pressure, P = — and it was proposed in [2] 
that the volume of the black hole be defined as the thermodynamic variable 
conjugate to P. Interpreting the mass as the enthalpy, equation (CQ) should 
then be modified to 

dM = TdS + VdP + VLdJ + $dQ, (2) 

where the thermodynamic volume is defined to be V = ^^\ SJ qi [TJ- The 
idea that A should be thought of as a thermodynamic variable that can be 
varied is not new and has been considered by a number of authors, [3]- [9]. 
Equation ([21), with Q = 0, was studied in [8], in the context of varying A, 
without any particular physical interpretation being given to the thermody- 
namic conjugate of A, often denoted G. 

One may question whether it is appropriate to identify A with a thermo- 
dynamic pressure. While a cosmological constant gives a pressure term in 
Einstein's equations, a fluid dynamical pressure is not necessarily the same as 
a thermodynamic pressure. In equilibrium situations however it is presum- 
ably correct to identify the fluid dynamical pressure with the thermodynamic 
pressure and we shall do so here. 

In general the enthalpy, H, is the heat energy beloved of chemists, it is not 
the internal energy, U, of the first law of thermodynamics. That distinction 
goes to the Legendre transform of the enthalpy, 

U = H-PV, (3) 
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where if is a function of S, P, J and Q while U (S, V, J, Q) is a function of 
purely extensive variables. Then we get the usual form of the first law, 



dU = TdS - PdV + VtdJ + §dQ. 



(4) 



This equation was written down in [9], using 6 — A notation, but its conse- 
quences were not pursued. 

When A is non-zero we should expect the PdV term to contribute to 
the mechanical energy that can be extracted from a black hole, by a Penrose 
process for example. For a negative A (positive pressure) the PdV term gives 
a positive contribution to dU if the black hole shrinks, and the PdV term 
reduces the amount of energy available for extraction as mechanical work W, 
with dW = —dU, hence reducing the efficiency. However we shall show that 
the maximal efficiency actually increases, relative to the A = case, when 
A < 0, because the maximal angular momentum of a black hole in AdS is 
greater than that of one with A = and this can outweigh the reduction 
in efficiency: the engine may not be as efficient at a given J but it can be 
pushed to higher J. Conversely one would expect that, for a positive A more 
energy becomes available at a given J, relative to A = 0, as the black hole 
shrinks. 

Of course there are no pistons pushing against a gas for a black hole, but 
a negative cosmological constant contributes a negative energy density to 
space-time so a shrinking black-hole exposes negative energy, thus increasing 
the black hole's internal energy and decreasing the amount of energy avail- 
able for mechanical work. Conversely a positive cosmological constant would 
presumably release extra energy as the black hole shrinks, that can be used 
to do work, hence increasing the efficiency at a given J. 

The most efficient way to extract energy from a black hole is in an isen- 
tropic process, with dS = and the area of the event horizon constant. 



We shall see that, for a rotating black hole, it is possible to reduce V while 
keeping S constant. The maximum amount of mechanical work that can be 
extracted in passing from an initial state i to a final state / is 



So 



dU > dU, 



mm 



-PdV + VtdJ + $dQ. 



(5) 




(6) 
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The efficiency is denned to be the ratio of the mechanical energy extracted 
to the initial heat energy (enthalpy), 



where the initial enthalpy is identified with the initial mass, Mj. It will be 
shown that, for A < 0, this can be as high as 52% for a rotating neutral black 
hole and 75% for a charged black hole. 

The volume of a black hole has only recently been considered as a thermo- 
dynamic variable, [2j [13] . At the simplest level, there is a natural tendency 
to assume that the area and the volume are related geometrically and are 
not independent. For a Schwarzschild black hole with radius r^, for exam- 
ple, the area is of course ^kt\ and indeed the thermodynamic volume works 
out to be Y r h but this seems co- incidental and no particular significance 
should be attached to it. At a deeper level it is not even clear how to define 
a volume geometrically, as the metric is not static in the interior of a black 
hole. For a Schwarzschild black hole which is not rotating we shall see that 
the thermodynamic volume and the area are not independent: fixing S fixes 
V so dV = in an isentropic process, and the PdV term does not contribute 
to the first law. But for a rotating black hole the area of the event horizon 
does not determine the thermodynamic volume uniquely and it is possible 
to vary the volume keeping the entropy constant, by changing the angular 
momentum and/or the charge. 

The properties of the thermodynamic volume and its contribution to the 
first law of black hole thermodynamics are explored in detail in this paper 
for a rotating charged black hole in four dimensional space-time with a neg- 
ative A. In section [2] thermodynamic potentials and the equation of state 
are discussed and the Legendre transform from the enthalpy to the internal 
energy is given explicitly. In section [3]the efficiency of a Penrose type process 
is analysed and section H] contains a discussion and outlook. Two appendices 
are dedicated to the technicalities of deriving some results used in the text 

2 The internal energy 

Including a pressure term in the first law gives 

dU = TdS - PdV + VLdJ + <S>dQ 
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where the internal energy, U(S, V, J, Q), is a function of extensive variables. 
The thermodynamic volume V is the conjugate variable to the pressure and 
is obtained from the mass, which is identified in |T] with the enthalpy, M = 
H(S,P,J,Q),by 

dH 



V dP 



(8) 

S,J,Q 



It was proposed in [2] that (jSJ) be defined to be the thermodynamic volume 
of the black hole. 

The line element for a charged rotating black hole in 4-dimensional anti-de 
Sitter space is [10] 



,9 A/, asm 2 9,\ p 2 , „ p 2 - A e sm 2 9 f , r 2 + a 2 , ,\ 

ds 2 = — - dt dcj) +*--dr 2 +^d6 2 +^- adt d(j) 

P 2 V - / A A e P 2 \ c J 

(9) 

where 

A = < r2 + ° 2 )f tfi -^W, A, = l-^cos^, 

a 2 

p 2 = r 2 + a 2 cos 2 ^, 5 = 1-—, (10) 

L l 

and the cosmological constant is A = — -p, which is related to the pressure 

by h - ¥■ 

The physical properties of this space-time are well known [11], and the 
first law, applied to this metric, was discussed in [12], but without a PdV 
term. The metric parameters m and g are related to the mass and charge by 

M =^> G = |- ( n ) 
The event horizon, r + , lies at the largest root of A(r) = 0, so 

u= <A±jWg£±i!L t (12) 

2r + L 2 

and the of the event horizon is 

A = iir (A±^l, (13) 
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The temperature is 



(L 2 + 3r^r 2 -a 2 ^ 2 -^)-^ 2 

4vrL 2 r + (r 2 +a 2 ) ' 1 ] 

The angular momentum, J = oM and the relevant thermodynamic an- 
gular velocity is 

L 2 (r+ + a 2 ) 



The electrostatic potential is 



r 2 +a 2 " 



One can scale L out from all the above expressions by defining dimen- 
sionless variables M = M/L, a = a/L, f+ = r + /L, etc., but we prefer to 
keep L explicit to expose more clearly the role of the pressure, and make the 
comparison with the L — > oo limit clear. 

Under the assumptions made here the thermodynamic volume for the 
Kerr-Newman-AdS black hole works out to be 



_ 2tt J (r 2 + a 2 )(2r 2 L 2 + a 2 L 2 - r 2 a 2 ) + L 2 q 2 a 2 
V = Y { L 2 E 2 r+ 

which is a simple generalisation of the Kerr-AdS volume derived in [13]. A 
direct derivation of ( 1T7|) from f lT2|) is most easily achieved by first writing the 
mass as a function of (S, P, J, Q), differentiating with respect to P, and then 
transforming back to (r+, a, q, L). 

When a = the area ( lT3"j) and volume (fT7|) are not independent and the 
area determines the volume uniquely, but when a is non-zero the area and 
the volume become independent. 

For asymptotically flat space, with L — > oo, one has 

v _ 2ir (r 2 + + a 2 )(2r 2 + + a 2 )+q 2 a 2 
3 r + 

While this reduces to the naive result, — for the Schwarzschild black hole, 
a geometrical interpretation when a is non-zero is not so clear. 

A correct description of the thermodynamics of the black hole, in terms of 
thermodynamic potentials, requires replacing the geometric variables (r + , L, a, q) 
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with thermodynamic variables (S,P, J,Q). The relevant expression for the 
mass, and hence the enthalpy, was derived in [9], 



1 l(S + nQ 2 
H(S,P,J,Q) :=-♦' 



8PS 2 



) 2 + 4 7r2 (1 + ^) J 2 



irS 



(19) 



This generalises the Christodoulou-Ruffini formula [2] for the mass of a 
rotating black hole in terms of its irreducible mass, M irr . The irreducible 
mass for a black hole with entropy S is the mass of a Schwarzschild black 
hole with the same entropy, M 2 rr = j-. 

In terms of thermodynamic variables the temperature is, [9], 



T 



OH 
~dS 



J,Q,P 



8nH 



ttQ 2 8PS 



and the thermodynamic volume is 

2 



V 



dH 



S,J,l 



3nH 



S(S + irQ J 



S 



8PS 



4vr 2 ^ 



J 
S 
(20) 



2vr 2 J 2 



(21) 



which is manifestly positive. 

The Legendre transform U = H—PV gives the thermal energy, a function 
of purely extensive variables. The transform is evaluated in an appendix to 
be 



U{S,V,J,Q) 



w 

Air 



Z7T j V 7T 



+ r 



(22) 



\J\ 



3V 
Air 



S 



7T 



1 

31 2 



SQ 2 
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Dimensional analysis implies that U is only a function of three independent 
variables, since U -t \U when S -t \ 2 S, V -t X 3 V, J -> A 2 J and Q -> XQ. 

One must be careful taking the J — > limit of these potentials. In this 
limit the enthalpy 



H(S,P,0,Q) 



1 



ttQ 2 8SP 



(23) 
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is linear in P and the Legendre transform is singular: 



dP 



(24) 



S,Q 3 W 



is independent of P and so cannot be inverted to obtain P(V). Conversely, 
when J is zero, the Legendre transform of (1231) is 

u^Ji( 1+ ^.), (25) 



2 V 7T V S . 

which is independent of the volume and is not equal to the J — > limit of 
(|22|) . unless a constraint, (f^-) 2 = (f) 3 , is imposed. Indeed (|25|) gives the 
wrong J = temperature (unless Q = P = 0, in which case U = H and 
T = = ^r- is the correct Hawking temperature for a Schwarzschild black 
hole). To get the correct temperature from ( 122]) in the J — > limit we must 
take the partial derivative with respect to S before setting J to zero, and 
take note of the fact that (f^-) 2 = (f) when J = 0. This constraint can 
be derived from (1221) directly by observing that U is not different iable with 

respect to J at J = unless ((f^ 1 ) 2 - (f ) 3 ) (</ 2 + vanishes there. It 
can also be seen directly when a = in ( Tl3i) and ( fl71) . 

To derive the relation between the pressure and the volume in general we 
first define v := and s := -. Then ( |22|) becomes 

t/(s, u, J,Q) = -^ {^(s 2 + sQ 2 + 2J 2 ) - |J|vV -s 3 ) (J 2 + sQ 2 )} (26) 
with temperature 

1 dU _ \ J\ {3J 2 (2v 2 - s 3 ) + sQ 2 (5v 2 - 2s 3 )} v ( s2 + ^Q 2 + 6J 



T 



vr ds 2n s 4 ^{v 2 - s 3 )(J 2 + sQ 2 ) 

and pressure 



2tts 4 



(27) 



_ _ £^ _ 3v\J\ J 2 + sQ 2 _ _3_ 

F ~ 47T ft, " 47TS 3 V v 2 - S 3 StTS 3 ' + s v + 2 J J . (28) 

In the J — > limit | J\ and aA> 2 — s 3 must vanish together for finite T and 
P. The equation of state, in the form of the relation between the pressure, 
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the temperature and the volume, can be obtained by eliminating s between 
(j2?p and (I28p . In the limit J — > the temperature is 



sttv^p + v 2 ' 3 - q 2 

T ^ = 4^ ' (29) 

which is the correct temperature as a function of P, Q and s, as can be seen 
from dHU, with a = 0, using (JT3J, (JT7J) and J = oM. When Q = fl2H} 
reduces to the equation of state in [2]. 
The small J correction to f )29|) is 

SttP^ + ^ - Q V/ 3 - 6 J 2 /_J*_ 

47ru 5/3 + ^9/3' ^9/3 J' ^ 

The angular momentum term dominates at small volumes for any J > 0. 
The first few terms in an expansion for the pressure are 

P = JL + ^ + (J1_ W (31) 

2v i/3 8vri; 2 / 3 8tu; 4 / 3 Anv 2 \ v 10 / 3 ' v 10 / 3 J ' 1 ; 

(the virial expansion for Q = is developed in more detail in the second 
appendix.) 

The P-V diagram is plotted to this order in figure 1 with Q = and 
J = 1 (for clarity the horizontal axis in figure 1 is v 1 ^ 3 rather than v). There 
is a critical point at T c w 0.0413, P c w 0.00280 and v t 1/3 w 3.08. This critical 
point represents a second order phase transition, similar in nature to that of 
a van der Waals gas — it is the same critical point that was found in [5] by 
fixing the pressure and varying J (the authors of [9] use P = = 0.119, 
corresponding to L = 1, and find a critical value of J at J c = 0.0236 giving 
^jr = 0.00282). The corrections in equation (|3T|) become large when w is 
very small, but at the critical point the error in the pressure is of order 
■ TT 5 7 3 ~ 10^ 5 , or about 1% of P c . When v 1 ^ 3 < 2 the corrections are of the 



r 1 



same order as P and the virial expansion breaks down. In the limit when 
v is small relative to J, a large J expansion shows that the leading order 
behaviour is P ~ J^-, which is the same low v behaviour as figure 1 but 
reduced by a factor 2. 

Thus by taking the J — > limit carefully the temperature can be kept 
finite when v 2 — > s 2 , but the non-analyticity in 022p is still present in the 
second derivative with respect to S, which diverges when J = and v 2 = 
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s 3 . This has the important consequence that the heat capacity at constant 
volume, Cy = Tj (§§)yjQ> tends to zero when J = 0, as observed in [2]. 
Cy can be non-zero for J ^ 0, and it is plotted in figure 2 for Q = as 
a function of S and J, with V — 1. For comparison the heat capacity at 
constant pressure, Cp = 77 (|^ ) PJ q, is plotted in figure 3 for Q = as a 
function of 5 and J, with P = 1. Cp vanishes when T = and diverges 
when |^ =0. 

A full stability analysis was carried out in [9] and figure 4 shows the phase 
diagram, when Q = 0, in terms of S/L 2 and J/L 2 . The top curve, labelled I, 
is the T = limit, the region above and to the left of this curve is unphysical 
because T < there. Cp in figure 3 vanishes on curve I. Curve II is the 
curve on which Q = 1/L, below which the black hole can be in equilibrium 
with radiation rotating at infinity [15]. On this curve asymptotic space- 
time is conformal to the 2 + 1 dimensional Einstein universe rotating at the 
speed of light, the region below and to the right of this curve is of relevance 
to conformal field theories. Curve III is the curve marking the edge of local 
stability, the determinant of the Hessian of the Gibbs free energy is divergent 
on this curve. The heat capacity, Cp, diverges on curve IV, which is blown 
up in the right hand picture (this curve was not included in the analysis 
of reference [9]). Curve V marks the Hawking-Page phase transition, below 
this curve the black hole is globally stable. The three curves marking the 
Hawking-Page transition, the edge of local stability and the Q = 1/L curve 
all have similar asymptotic forms for large S/L 2 , -p ps 2 ^ 2 l a + ^ (z^) - 

In the next section it will be more convenient to use the variables S, P, 
J and Q, rather than S, V , J and Q, in terms of which 

(5 + vrQ 2 ) (5 + vrQ 2 + «p) + 47T 2 (1 + J 2 
U = - _ — , (32) 

2 A hs 



(S + nQ 2 +*rf-y + 4n 2 (l + ^) J 2 

which is manifestly positive. In terms of geometrical variables, 

_ (r 2 + a 2 )(2L 4 + a 2 r 2 - a 2 L 2 ) + g 2 L 2 (2L 2 - a 2 ) 
U ~ 4i> S 2 r+ • (33) 
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3 Extracting energy from rotating black holes 



For an electrically neutral black hole, we set Q = in ( 1321) and the internal 
energy U = H — PV, written in terms of S, L 2 = -^p and J, is 

rr /? 2 (1 + ^ +^(1 + A) ^ 



In an isentropic, isobaric process the black hole can yield mechanical work 
by decreasing the angular momentum. If J is reduced from some finite value 
to zero the efficiency is 

_ U(J)-U(0) 



H(J) 

^(1 + ^+4^(1 + ^) J 2 

S 



(35) 



For a given S and L there is a maximal value of J determined by de- 
manding that T > 0, which requires 

^<S( 1 + ^)( 1 + §> (36) 

The greatest efficiency is for extremal black holes, when the bound (13"6"j) is 
saturated, 

' 2(27rL 2 + 3,S) (nL 2 + S) V2nL 2 + 3S' 1 ' 

The asymptotically flat case (A = 0) is obtained by taking the limit L — > oo, 
77 = 1 — 4= = 0.29 . . ., which is the familiar result [16]. For finite L the effi- 
ciency is greater than the asymptotically flat value. It has a maximum value 
of 0.5184... when S/L 2 = 15.39... (obtained by solving a quartic equation for 
S), and asymptotes to 0.5 for large S. 

Differentiating ffl9|) with respect to P, we see that the thermodynamic 
volume is a monotonically increasing function of J, hence the volume de- 
creases as J is lowered, keeping the area constant. The black hole also heats 
up during the process. 
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When the black hole is charged there are various possibilities for ex- 
tracting work through an isentropic, isobaric process: one could decrease J 
keeping Q constant or decrease Q keeping J constant or decrease both si- 
multaneously. Since all that matters is the initial and final values of U, the 
initial and final values of J and Q uniquely determine the efficiency in such 
a process. 

If both J and Q are both non-zero initially and are both decreased to 
zero, the efficiency is 

(S + vrQ 2 ) (S + vrQ 2 + + 4vr 2 (l + ^) J 2 

jj — i L (28) 

(S + vrQ 2 + ^) 2 + 4^(1 + ^) J 2 

S 



(S + *Q 2 + S) + 4vr 2 (1 + ^) J 2 

The optimal efficiency is for extremal black holes and the extremal value of 
J is given by 



J 2 =i 

max 



\4 



(39) 



allowed range of Q 2 



when T in ( 1201) vanishes. Positivity of J max then imposes a restriction on the 

")2 



°^(f)( 1 + §)- < 40) 

The greatest efficiency is for an extremal black hole with maximal charge, 

^ © + §) • < 41) 

when 

V = nL2 + 3S (42) 

V 2(vrL 2 + 2S)- 

For large S efficiencies of up to 75% are possible in principle for extremal 
charged, rotating AdS black holes, though as we saw in the previous section 
this is reduced to 51.8% for extremal rotating, but electrically neutral, black 
holes. For comparison, the equivalent figures in asymptotically flat space- 
time, with A = 0, are 50% and 29.3% respectively. 
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4 Conclusions 



The thermodynamics of rotating black holes in four dimensional in anti- 
de Sitter space-time has been discussed in detail, with particular attention 
payed to the role of pressure and the volume. The negative cosmological con- 
stant is interpreted as a positive pressure and treated as a thermodynamic 
variable whose conjugate thermodynamic variable is a thermodynamic vol- 
ume. For rotating black holes the thermodynamic volume is independent of 
the geometric area, and hence independent of the entropy, allowing a full 
thermodynamic treatment. 

The black hole mass is associated with the enthalpy H, as suggested in 
PP , and the internal energy U = H — PV is lowered relative to the mass by 
a positive pressure. Vanishing angular momentum is a non-analytic point of 
the thermodynamics and must be treated with some care, in particular one 
cannot fully understand black hole thermodynamics simply by focusing on 
the non-rotating, Schwarzschild case — rotation is essential for a complete 
analysis. 

The black hole equation of state has been analysed in terms of pressure 
and volume and a virial expansion developed. Non-zero angular momentum 
causes a rapid rise in pressure as the volume is reduced at constant tempera- 
ture, a rise that is not present at zero angular momentum. The critical point 
found in p| is found to be very similar in form to that of a van der Waals 
gas. 

In a Penrose process the volume decreases as the black hole looses angular 
momentum and the PdV term in the first law reduces the amount of energy 
available to do work. However a negative cosmological constant also extends 
the range of the angular momentum, pushing the maximum allowed value 
beyond that of the A = case, so that in fact more work can be extracted 
from an extremal black hole in a space-time that is asymptotically anti-de 
Sitter than from one in a space-time that is asymptotically flat. Efficiencies 
of up to 75% are theoretically achievable for a charged black hole in asymp- 
totically anti-de Sitter space-time, compared to 50% in asymptotically flat 
space-time. For electrically neutral black holes the corresponding figures are 
51.8% compared to 29.3%. 

The analysis here has focused on the A < case for three reasons. Firstly 
there is an ambiguity in the definition of thermodynamic variables in the 
A > case, there are two relevant event horizons to be considered and hence 
two different temperatures in general. For a general choice of black hole 
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parameters it is not possible to make the Euclidean geometry regular by a 
unique periodic identification of Euclidean time. Secondly A > implies 
negative pressure and consequent thermodynamic instabilities. This latter 
problem is not necessarily a fatal objection, black holes are well known to be 
unstable for A = because the heat capacity, more precisely Cp, is negative, 
but we can still consider metastable situations in which a temperature can 
be defined for a period of time much shorter than the timescale of thermal 
instability. Indeed there are regions of parameter space for A < in which 
Cp < 0, and we can still learn something about physical properties of black 
holes in this regime. Negative pressures can be useful in metastable situations 
[T?] . Lastly it is not even clear how to define the mass of a black hole in a 
A > space-time as there is no accessible asymptotic regime in which one 
can compare with the zero mass solution. 



Appendix 1 

To calculate the internal energy U = H — PV, we write the enthalpy in the 

form 

H= ^a + pP + 7P 2 , 



where 

a : = 

6 := 

7 := 

Note that the discriminant, 



j 



Air (S fS 



3 [ IX \7T 



Q 2 + 2f 



/3 2 - 4a~f 



J 2 J 2 



71 



is positive. 
Now 



V 



dH 



S,J,Q 



(3 + 2-fP 
2H 



P 



2HV-/3 
27 



(43) 



(44) 



(45) 



(46) 
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This allows us to re-express H function of V, 



H=\r- ia \ (47) 



We can immediately conclude that 



V > ( f ) (f '] , (48) 



in agreement with the observation in [15] . 
It is now straightforward to determine 



U = H-PV = H-[^ _gLg_vggEM, (49 ) 



7 27 / 27 27 
which immediately gives (122]) in the text 



Appendix 2 



In this appendix the virial expansion is developed. For simplicity we set 
Q = 0, but the same techniques can be applied to the case of non-zero Q. 
To develop the expansion we use dimensionless variables 

v s 8ttPJ 



t := 2nTJ 1/2 , (50) 



J 3 / 2 ' J' 3 

in terms of which equations ( )2Tj) and ( |28]) can be written 

^ = 3(2y 2 -x 3 ) _ y(x 2 + 6) 

2y x 2 + 2 
P = r~F= = — ■ (52) 

When y and x are large let y 2 ~ x 3 , with y 2 — x 3 = z 2 , then t is finite 
provided 2 ~ in which case p ~ ^73- Replacing a; with 2; in equations 
(IBT]) and gives 

3(y 2 + z 2 ) = {t(y 2 -z 2 Y /3 + y(y 2 -z 2 f 3 + 6y}z, (53) 

(y 2 -z 2 )zp = 2y- z\(y 2 -z 2 ) 2 ^ + 2 



2 1 

P ~ ^(|/ + z) ~ (y 2 -z 2 y/ 3 ' ^ ' 
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Now we expand in powers of u = -373 . Let z = ^rC for some ((t, u), in terms 
of which (ED) becomes 



t 2 / I V t 2 / V t 2 

(55) 

from which we can immediately conclude that 

C= , * 5 +u(t,u), (56) 

where an expansion of u in u starts at order 10. Any desired order can be 
obtained by further expanding u 

00 

,10 



5>„(tK, (57) 



00 = u 

n=0 

with the co-efficients a n (t) to be determined. Putting (1571) and (1551) into 
(155j) . and equating co-efficients of powers of tt, gives an iterative procedure 
for evaluating the co-efficients a n which can then be used to show that: 

u 10 it 11 u 12 u 13 u u 2 + 3*4 lfi , , . 

w = V ~ 1 V + 2 V " 4 V + 5 V 36 + ( u ) • (58) 

Finally using this expansion in ( )54|) gives the required virial expansion 

2tu 1 9 K in / it u 2 u 3 u 4 u 5 . , 

p = — -5" + 2 « -« ( 8 7- 9 F + 1 V _1 V + 1 V (59) 



. 13 i±^„«) +0 K') 



the first three terms of which are used in the text. 
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Figure 1: P-V diagram for J = 1 and Q = 0, plotted using the approximation 
in equation (13TT) . P is plotted as a function of vs for T = 0.02, 0.025, 0.3, 
0.35, 0.04, 0.045 and 0.05. The critical point for this value of J is T c w 0.0413, 
P c ~ 0.00280 and f c ~ 3.08. J > causes the rapid rise in P at low values of v. 
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Figure 2: Cy as a function of S/L 2 and J/L 2 , with fixed volume set to L 
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Figure 3: Cp as a function of S/L 2 and J/L 2 , with fixed pressure 
set to J^, i.e L = 1. Cp diverges along curve IV in figure 4 and 
vanishes along curve I. 
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Figure 4: phase diagram for q = 0, plotted in terms of S/L 2 and J/L 2 . 
The region above I is forbidden, because T < 0; 

in the region above curve II the 3-d Einstein universe at infinity rotates 
faster than the speed of light; 

curve III bounds the region of local stability (an analysis of the Gibbs free 
energy shows that the black hole is locally unstable above curve III); 
Cp diverges on curve IV; 

in the region above curve V the black hole is unstable due to the Hawking- 
Page phase transition; 

in region below curve V the black hole is stable. 
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